UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADPO 14047 

TITLE:  Stochastic-Constraints  Method  in  Nonstationary  Hot-Clutter 
Cancellation  Part  II:  Unsupervised  Training  Applications 

DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 
Availability:  Hard  copy  only. 


This  paper  is  part  of  the  following  report: 

TITLE:  Military  Application  of  Space-Time  Adaptive  Processing  [Les 
applications  militaires  du  traitement  adaptatif  espace-temps] 

To  order  the  complete  compilation  report,  use:  ADA415645 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  etc.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 

ADPO  14040  thru  ADPO  14047 


UNCLASSIFIED 


8-1 


Stochastic-Constraints  Method  in  Nonstationary  Hot-Clutter 
Cancellation  — Part  II:  Unsupervised  Training  Applications 

Prof.  Yuri  I.  Abramovich 1,2 

1 Surveillance  Systems  Division,  Defence  Science  and  Technology  Organisation  (DSTO) 
PQ  Box  1500,  Edinburgh  S.A.  5111,  Australia 
2 Cooperative  Research  Centre  for  Sensor  Signal  and  Information  Processing  (CSSIP) 
SPRI  Building,  Technology  Park  Adelaide,  Mawson  Lakes  S.A.  5095,  Australia 

yuri@cssip.edu.au 


Abstract 

We  consider  the  use  of  “stochastically  constrained”  spatial  and  spatio-temporal  adaptive  processing  in 
multimode  nonstationary  interference  (“hot  clutter”)  mitigation  for  scenarios  that  do  not  allow  access  to  a 
group  of  range  cells  that  are  free  from  the  backscattered  sea/terrain  signal  (“cold  clutter”).  Since  super- 
vised training  methods  for  interference  covariance  matrix  estimation  using  the  cold-clutter-free  ranges  are 
inappropriate  in  this  case,  we  introduce  and  analyze  adaptive  routines  which  can  operate  on  range  cells 
containing  a mixture  of  hot  and  cold  clutter  and  possible  targets  (unsupervised  training  samples).  Theoret- 
ical and  simulation  results  are  complemented  by  surface-wave  over-the-horizon  data  processing,  recently 
collected  during  experimental  trials  in  northern  Australia. 

1.  Introduction 

In  this  lecture,  we  continue  our  study  into  the  adaptive  spatial  and  spatio-temporal  mitigation  of  nonstationary 
interference  in  the  field  of  high  frequency  (HF)  over-the-horizon  radar  (OTHR).  The  preceding  lecture  dis- 
cussed the  fundamental  ideas  behind  the  new  “stochastic-constraints”  approach,  which  has  been  proposed  to 
achieve  effective  hot-clutter  suppression  whilst  maintaining  distortionless  output  cold-clutter  post-processing 
stationarity.  The  introduced  operational  implementation  is  based  on  the  availability  of  range  cells  specifically 
free  of  cold  clutter,  a situation  that  is  quite  typical  of  pulse-waveform  (PW)  OTHR. 

Meanwhile,  even  within  the  PW  OTHR  architecture,  attempts  to  increase  the  radar  duty  cycle  (ie.  the  trans- 
mitted energy)  usually  lead  to  multifrequency  operation  with  a transmission  pause  for  each  given  frequency, 
limited  to  the  operational  range  depth.  For  surface-wave  (S  W)  OTHR,  such  a multifrequency  mode  of  operation 
leaves  little  room  for  (sea-)  clutter-free  ranges  being  available.  For  frequency-modulated  continuous  waveform 
(FMCW)  OTHR,  the  accessibility  of  such  range  cells,  containing  interference/jamming  and  noise  signals  only, 
is  also  problematic.  Traditionally,  FMCW  systems  operate  with  a linear  frequency-modulated  (LFM)  waveform 
and  have  receiving  systems  which  employ  a mixing  of  the  received  waveform  with  a time-delayed  version  of  the 
transmitted  waveform.  This  process  is  known  as  “deramping”  and  is  followed  by  spectral  analysis  (weighted 
Fourier  transform)  within  the  comparatively  narrow  bandwidth  of  a low-pass  filter  at  the  output  of  the  mixer. 
This  bandwidth  is  usually  adjusted  to  the  operational  range  depth,  so  that  the  frequencies  corresponding  to 
the  “skip  zone”  ranges  are  usually  filtered  out  at  this  stage.  Thus  only  “operational”  range  cells  which  con- 
tain interference,  cold  clutter  and  possible  targets  (unsupervised  training  samples)  are  available  for  any  type  of 
processing. 

Note  that  the  “simple”  solution  of  increasing  the  low-pass  filter  bandwidth  has  its  limit  because  of  the 
significant  increase  in  range  sidelobes,  due  to  the  frequency  mismatch  [1].  Unfortunately,  these  range  sidelobes 
also  smear  the  multimodal  interference  structure  somewhat. 

Obviously  the  type  of  operational  routines  discussed  in  [2]  based  on  access  to  cold-clutter-free  ranges  are 
not  applicable  to  these  types  of  OTHR.  It  is  important  to  emphasize  that  most  existing  adaptive  algorithms  are 
unable  to  address  the  problem  of  effective  interference  and  cold-clutter  mitigation  simultaneously.  Indeed,  the 
number  of  “sweeps”  (or  repetition  periods,  N ) within  the  coherent  processing  interval  (CPI,  “dwell”),  combined 
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with  the  number  of  spatial  receiving  channels  M,  is  typically  too  large  to  consider  even  slow-time-only  STAP 
as  a solution  to  this  problem,  since  the  problem  has  MN  degrees  of  freedom,  and  we  have  M ~ 32  to  64  and 
N ~ 128  to  256  for  typical  radar  facilities.  Not  only  is.  this  solution  impractical  for  computational  reasons,  but 
the  number  of  available  range  cells  T within  the  bandwidth  of  the  above-mentioned  low-pass  filter  is  usually 
too  small  ( T — 50  to  70)  compared  with  the  product  M N to  provide  reliable  convergence  for  the  STAP  filter 
to  the  efficiency  given  by  the  true  (deterministic)  covariance  matrix  as  T — > oo. 

Moreover,  this  problem  is  aggravated  if  severe  multimode  ionospheric  propagation  [3,  4,  5,  6]  gives  rise  to 
the  so-called  “hot-clutter”  phenomenon,  whereby  a signal  produced  by  a single  interference  source  is  seen  as  a 
multiplicity  of  interference  signals  at  the  receiving  array.  According  to  [4]: 


Ionospherically  propagated  signals  may  consist  of  several  modes,  nonstationary  in  bearing  and 
highly  correlated  with  each  other.  Furthermore,  due  to  the  signal  being  refracted  from  an  inhomo- 
geneous region,  each  mode  can  be  considered  to  consist  of  a specular  ray  surrounded  by  a cone  of 
diffracted  rays.  The  resulting  wavefront  for  each  mode  may  therefore  be  far  from  pi  an  sir. 


When  the  temporal  support  of  such  multimode  is  significant,  fast-time  STAP  with  Q consecutive  range  bins 
(“taps”,  Nyquist  intervals)  involved  in  adaptive  processing  is  one  way  to  enhance  hot-clutter  rejectability.  In 
some  situations  this  is  the  only  way  to  reject  multimode  interference,  especially  if  one  of  these  modes  affects  the 
main  beam  of  the  receiving  array  [2],  Obviously  if  MN-v  ariate  slow-time  STAP  is  considered  to  be  impractical, 
this  is  much  more  the  case  for  the  “fully  adaptive”  MA^Q-variate  STAP  described  in  [7,  8]. 

Thus  from  a practical  viewpoint,  we  should  consider  a scheme  whereby  each  “finger  beam”  is  associated 
with  an  M -variate  spatial-only  adaptive  process  (SAP),  or  at  most  with  an  MQ -variate  fast-time  STAP  for 
interference-only  (hot-clutter)  mitigation.  The  output  signals  of  each  beam  should  then  be  processed  by  some 
standard  slow-time  inter-sweep  coherent  processing  technique.  Most  existing  HF  OTHR  facilities  adopt  a 
weighted  Fourier  transform  for  the  detection  of  moving  targets  in  the  Doppler  frequency  domain.  Clearly 
these  standard  cold-clutter  post-processing  techniques  can  only  be  effective  if  the  properties  of  the  (scalar) 
cold  clutter  signal  at  the  output  of  the  adaptive  beamformer  are  not  contaminated  by  beamformer  fluctuations 
caused  by  adaptation.  This  consideration  is  of  great  importance,  since  spatial  and  spatio-temporal  properties  of 
the  ionospherically  propagated  jamming  signal  (hot  clutter)  are  highly  nonstationary  over  the  CPI.  One  possible 
analytic  model  that  involves  variations  of  the  ionospheric  propagation  medium  during  the  CPI  is  introduced  in 
[6],  and  further  mentioned  in  [2]. 

While  the  interference  mitigation  efficiency  of  an  adaptive  filter  that  is  constant  over  a CPI  may  be  quite 
poor  in  such  cases,  we  may  achieve  a substantial  improvement  only  if  the  corresponding  SAP  or  STAP  system 
can  track  perturbations  of  the  hot-clutter  covariance  matrix  along  the  CPI.  When  perturbations  of  the  antenna 
array  pattern  due  to  the  tracking  are  unconstrained,  we  observe  a significant  degradation  in  the  properties  of 
the  cold-clutter  Doppler  spectrum  (peak-to-sidelobe  ratio,  “subclutter  visibility”)  [2].  It  has  been  demonstrated 
[9,  10,  11,  12,  2]  that  additional  data-dependent  (“stochastic”)  constraints  could  be  imposed  on  adaptive  inter- 
ference rejection  filters  to  maintain  distortionless  output  cold-clutter  processing  stationarity. 

Here  we  introduce  and  analyze  operational  stochastically  constrained  (SC)  SAP/STAP  routines  that  im- 
plement the  same  fundamental  principle,  but  for  unsupervised  training.  The  basic  idea  of  the  first  introduced 
approach  is  straight-forward:  we  reject  the  cold  clutter  using  a moving-target  indicator  (MTI),  calculate  the 
STAP  weights  from  this  data  that  will  reject  hot  clutter  but  maintain  cold-clutter  stationarity,  then  apply  these 
weights  to  the  original  data,  enabling  the  cancellation  of  cold  clutter  using  standard  Doppler  processing. 

When  the  MTI  filter  can  be  designed  in  a nonadaptive  fashion,  the  implementation  of  this  basic  idea  is 
quite  simple,  and  appears  to  be  already  discussed  in  classified  papers  [13,  14,  15].  Therefore  we  concentrate 
on  the  efficiency  analysis  of  this  approach  for  typical  HF  OTHR  scenarios.  In  the  more  general  case  when  the 
cold-clutter  spectral  properties  are  unknown  a priori,  implementation  of  the  basic  idea  is  less  trivial,  and  the 
hot-clutter-only  rejection  filter  should  be  properly  retrieved  from  the  associated  3-D  STAP. 
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2.  Operational  Routine  for  Unsupervised  Training 

Our  goal  is  now  to  construct  a complementary  operational  routine  which  implements  an  approximated  version 
of  the  generic  solution  for  the  unsupervised  training  scenario.  We  may  formulate  a similar  approach,  providing 
we  can  somehow  obtain  proper  rejection  of  the  hot-clutter  component.  Clearly  we  must  be  prepared  to  sacrifice 
some  performance  compared  with  supervised  training,  due  to  the  absence  of  samples  containing  only  hot  clutter 
and  noise. 

2..1  Known  (Local)  Cold-Clutter  Model 

In  [2],  we  demonstrated  using  real  OTHR  sea-clutter  data  that,  despite  the  quite  complicated  nature  of  the  sea- 
clutter  Doppler  spectrum  over  the  entire  dwell,  stabilization  by  stochastic  constraints  can  be  quite  effective, 
even  if  we  adopt  a simplified  local  low-order  AR  model  (ie.  accurate  only  for  a limited  number  of  consecutive 
sweeps).  In  fact,  it  has  long  been  known  that  a simplified  one-stage  MTI  filter,  for  example,  can  effectively 
reject  cold  clutter  even  if  this  clutter  is  much  more  complicated  than  an  ordinary  Markov  chain,  when  such  a 
filter  is  truly  optimal. 

According  to  the  general  philosophy  discussed  in  Section  I,  our  goal  is  to  apply  the  “shortest”  possible  MTI 
filter  to  reject  the  backscattered  cold-clutter  first,  in  order  to  extract  the  hot-clutter  samples  that  can  be  later  used 
for  hot-clutter  covariance  matrix  sample  averaging.  Here  it  is  crucial  to  involve  the  minimum  possible  number 
of  repetition  periods,  due  to  the  nonstationarity  of  the  hot  clutter,  which  increases  the  signal  subspace  dimension 
of  the  averaged  covariance  matrix  R^+1.  On  the  other  hand,  if  the  cold-clutter  residues  are  significant,  the  hot- 
clutter  covariance  matrix  estimate  will  be  significantly  compromised. 

Since  the  AR  model  for  SC  purposes  also  needs  to  be  of  minimum  possible  order  (/c),  it  could  be  used 
directly  to  design  the  MTI  filter  providing  that  the  parameters  bj  (j  = 1, . . . , k)  and  <r|  are  known  a priori  or 
estimated.  Thus  we  assume  that  a (k+  l)-variate  slow-time  preprocessing  filter 

z%t  = Zkt  + y>;  Zk-j,t  for  k = K+l,...,N  (1) 

j- 1 

would  lead  to  effective  cold-clutter  mitigation  within  the  output  snapshot  zkt> ie- 

£ {~zlt  *$}  - <?l  + X)  M2  &b2j  + <4  Imq  . (2) 

j-0 

with  b0  — 1,  and 

£lM2*K  • (3) 

j— o 

This  gives  rise  to  the  following  straight-forward  operational  routine. 


Algorithm  1 

Step  1 Given  the  observed  M-variate  snapshots  Zkt  (k  — 1, . . AT;  t = 1 , . . T),  form  the  “stacked” 
AfQ-variate  samples  Zkt  (k  — 1, . . N;  t — 1, . . . ,T—Q+1). 

Step  2 Preprocess  the  samples  z&  by  the  (k+1) -variate  slow-time  preprocessing  filter  with  impulse  response 
[1,  bi, . . . , bK]  to  get  the  cold-clutter-free  training  samples  z£t  using  (1). 


Ro  < 


Step  3 Compute  the  sample  estimate  of  the  hot-clutter-plus-noise  averaged  covariance  matrix 


Rk 


1 

T-Q+l 


T—Q+l 

E~k  ~kH 
zkt  zkt 


t= i 


for  k — /c+1,  • . N ■ 


(4) 


Step  4 Construct  the  initial  (k  = k + 1)  STAP  filter 

<V+I,t  = (Cl)"1  A^o)  [if(^o)  (Ci)_1^w]  (5) 

Note  that  this  filter  is  stochastically  unconstrained  and  so  is  range-invariant  x). 

Step  5 For  the  next  adaptive  filter  (now  range-dependent  and  stochastically  constrained),  apply  “slid- 

ing window”  averaging  to  the  preprocessed  input  data  z*+2/. 


~ 0.1/ 

Rk+2 


1 

T-Q+l 


T-Q+l 
t= 1 


-K.H 


zZ 


k+2  ,t 


in  order  to  ensure  hot-clutter  rejection  over  the  repetition  periods  k = 2,...,  k+2. 
Compute  all  subsequent  STAP  filters  (k  = k+2 , . . . , N)  by 


~ av 
wkt 


(6) 


(7) 


where 

ZQV 

Rk  = 

and  the  system  of  stochastic  constraints  is 


1 

T-Q+l 


T-Q+l 

£ 


t= i 


where 


Aq(0 0) 


*(go) 

o)J 


Zkt 


— \zk~l,t 


zk-2,t 


Zk-K,t\ 


(8) 


(9) 

(10) 


Note  that  the  particular  parameters  used-in  [9,  2]  to  simulate  HF  scattering  from  the  sea  - 

k=2,  b0  = l,  bi  = -1-9359,  b2  = 0.998 , <r^  = 0.009675  (11) 

means  that  all  three  covariance  matrices  R^  (j  = 1,2,3)  are  almost  equally  weighted  by  \bj\2,  while  the 
cold-clutter  component  in  zk+11  is  rejected  to  the  level  of  the  weak  temporally  white  noise. 

a ax' 

In  Step  5,  note  that  the  k repetition  periods  k = 2, ...,«+ 1 are  common  in  the  construction  of  7?,.+1  and 

to V 

Rk+ 2-  The  system  of  k stochastic  constraints  corresponding  to 

™ktzk-j,t  = wk_uzk-j,t  for  j (12) 


may  then  be  written  as 


™K+i,t  zkt  = wa:+2,t  zkt  for  k=  2, ...,«+ 1 . 


(13) 
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It  is  worth  emphasizing  that  the  hot-clutter  covariance  matrix  estimates  are  computed  using  the  preprocessed 
samples  zfct,  while  the  stochastic  constraints  use  the  original  snapshots  z & . Thus  the  right-hand  sides  of  these 
constraints  consist  of  the  cold-clutter  samples  mainly  due  to  the  properties  of  the  estimate  (6),  while  condition 
(3)  ensures  that  the  filter  wf+2j  properly  processes  the  cold-clutter  samples,  since  the  hot  clutter  is  to  be 
rejected. 

Our  final  comment  is  a reminder  that  the  matrix  in  the  square  brackets  in  (5)  has  dimension  q X q,  while  that 
of  (7)  is  (q+n)  x (q+n),  since  all  STAP  filters  but  the  first  have  k additional  data-dependent  linear  constraints. 

A comparison  of  this  algorithm  with  its  counterpart  for  supervised  training  (the  subject  of  the  previous 
lecture),  reveals  that  the  only  significant  difference  between  these  two  techniques  is  how  we  obtain  the  cold- 
clutter-free  training  samples  to  estimate  the  hot-clutter-plus-noise  covariance  matrix.  For  unsupervised  training, 
we  heavily  rely  on  the  ability  of  the  K-order  slow-time  filter  rejecting  the  cold-clutter  signal  practically  to  the 
level  of  white  noise.  In  this  regard,  we  stress  that  such  rejectability  validates  the  use  of  a low-order  AR  cold- 
clutter  model.  On  the  other  hand,  suitability  of  a low-order  AR  model  to  real  cold  clutter  over  a very  limited 
number  of  repetition  periods  (ie.  locally)  does  not  mean  that  the  same  model  is  appropriate  for  describing  cold 
clutter  over  the  entire  CPI  (ie.  globally).  Therefore  the  proposed  preprocessing  cannot  be  used  as  a moving 
target  indicator  (MTI)  filter  as  such  since,  for  low-speed  targets,  this  processing  is  far  from  being  globally 
optimal  and  leads  to  an  unacceptably  broad  range  of  “blind”  Doppler  frequencies.  It  is  thus  essential  that 
the  proposed  preprocessing  is  applied  only  to  extract  hot  clutter  which  is  uniformly  distributed  over  the  entire 
Doppler  frequency  band,  while,  for  proper  target  detection,  we  retain  the  output  scalar  cold-clutter  signal  as 
close  as  possible  to  its  original  form. 

It  was  demonstrated  in  [2],  by  processing  of  real  skywave  and  surface-wave  OTHR  data,  that  the  low- 
order  AR  model  supports  the  stochastic-constraints  approach  by  preserving  the  global  structure  of  the  output 
scalar  cold-clutter  signal,  which  is  more  complicated  than  suggested  by  the  model.  In  Section  V,  we  shall 
demonstrate  that  the  low-order  AR  model  also  provides  effective  real  cold-clutter  rejection  by  the  associated 
preprocessing  filter.  Nevertheless,  none  of  the  preprocessing  schemes  can  entirely  reject  real  cold  clutter  and  so 
some  cold-clutter  residues  are  always  present  within  the  training  sample  . This  means  that  additional  losses 
are  expected  compared  with  supervised  training,  especially  if  these  residues  are  above  the  white  noise  level. 

Another  source  of  loss  stems  from  a reduction  in  the  training  sample  size:  for  supervised  training  we 

,XT) 

separately  estimate  each  Rk  using  the  set  of  hot-clutter  training  samples  within  each  repetition  period  k, 
while,  for  unsupervised  training,  we  use  only  one  set  of  (T — Q+ 1)  ranges  to  estimate  a single  matrix  which 
approximates  the  sum  of  all  («+ 1)  matrices  . 

2..1  Known  Order  of  the  (Local)  Cold-Clutter  Model 

In  some  applications  it  is  reasonable  to  admit  that  a preprocessing  MTI-type  filter  can  be  defined  based  on 
some  fairly  standard  a priori  assumptions.  In  the  field  of  SW  OTHR,  for  example,  such  a local  model  might 
be  inferred  from  an  analytical  description  of  the  sea-echo  Doppler  spectra  [16].  In  general,  though,  our  prior 
information  is  much  more  limited.  To  accommodate  this,  we  specify  only  the  AR  model  order  k prior  to 
data  processing,  while  the  actual  parameters  values  bj  are  unknown.  Therefore,  both  the  cold-clutter  temporal 
properties,  as  well  as  the  hot-clutter  spatial  (for  SAP)  or  fast-time  spatio-temporal  (for  STAP)  properties,  are  to 
be  estimated  simultaneously  using  the  unsupervised  data  zkt.  It  should  be  clear  that  this  can  be  achieved  only 
by  incorporating  slow-time  to  deal  with  the  cold  clutter. 

The  main  idea  behind  the  algorithm  that  we  propose  is  again  straight-forward.  For  the  properly  chosen 
AR  model  of  order  n,  the  appropriate  MQ(k+  1) -variate  3-D  STAP  filter  wkt  should  then  effectively  reject 
both  hot  and  cold  clutter'  On  the  other  hand,  since  the  hot-clutter  signal  is  uncorrelated  over  adjacent  sweeps, 
effective  mitigation  will  occur  only  if  each  of  the  («+l)  fast-time  MQ-variate  2-D  STAP  subfilters  effectively 
rejects  hot  clutter.  Thus,  using  the  first  MQ  components  of  the  3-D  STAP  filter,  we  may  reject  the  hot-clutter 
component,  while  the  cold-clutter  contribution  should  remain  unperturbed.  Clearly  we  are  decomposing  the 
overall  3-D  STAP  filter,  extracting  the  component  responsible  for  hot -clutter  mitigation  only. 
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Let  us  introduce  the  MQ(k-\-  1) -variate  “doubly  stacked”  vector 


*kt  = 


Zk—K,t 

Zk-l,t 

Zkt 


In  the  absence  of  a target  signal,  we  may  write 

£{'nk}  = ^ {*****!}  = + £{yktVkt}  + ^Jmq(k+ i) 

since  the  hot-  and  cold-clutter  signals  are  mutually  independent.  Correspondingly, 


| = s{xktxkt}  + ^Imq{k+1)  = 


R 


XT] 

k—K 


R 


XT) 

fc-K+1 


nr 


(14) 


(15) 


(16) 


since  the  hot  clutter  is  uncorrelated  over  adjacent  repetition  periods.  For  a cold-clutter  component  described  by 
the  scalar  multivariate  model 

07) 

j=x 

we  have 

j = e{htn<}  = Kk+i®#?  (is) 

where 

nK+i  = Toep[r0,...,r«]  (19) 

rj  (j  — 0, . . . , k)  are  the  inter-period  (scalar)  correlation  lags,  and  Ry  is  defined  by  the  correlation  property 

£ { Vkt  Vk't ' } = S(t.  - f)  Ryk_k,  (20) 

and  by  (2). 

In  order  to  justify  some  of  our  assertions,  let  us  temporarily  simplify  our  problem.  Suppose  we  ignore  the 
nonstationarity  of  the  hot  clutter  over  the  interval  of  (k+1)  adjacent  repetition  periods,  and  let  us  assume  that 
the  cold-clutter  spatial  covariance  matrix  is  diagonal, le.  Ry  = Imq  . (The  angular  distribution  of  cold  clutter  is 
always  assumed  to  be  sufficiently  wide  to  exclude  any  possibility  that  pure  SAP  can  solve  the  problem;  ftirther 
discussion  appears  in  [2].)  Thus  the  covariance  matrix  Ry  is  always  well-conditioned,  and  so  this  second 
assumption  does  not  change  the  essence  of  the  problem. 

Under  these  two  simplifying  assumptions,  we  may  present  the  overall  MQ(k-\-\)  -variate  covariance  matrix 
- z 

Rk  in  the  form 

^{zktZkt}  = ttk  = R+i^RT  + UK+1  <8>  Imq  ■ (21) 

It  is  known  [17]  that  if  Xf  and  uf  (s  = 1, . . m)  are  respectively  the  eigenvalues  and  eigenvectors  of  some  m- 
variate  matrix  A,  while  and  fiB  (t  = 1 , . . . , n)  correspond  to  some  n -variate  matrix  B,  then  the  eigenvectors 
uf  (j  = 1, mn)  of  the  matrix 

C = Im  <g>  B + A ® In  (22) 

are  given  by 

uf  — uf®uf  for  s = 1, . . .,m;  t = 1, . . ,,n  (23) 
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and  correspond  to  the  eigenvalues  of  C equal  to  (Xf  -f  Af ). 


By  this  property,  the  MQ(/c+l) -variate  vector  wk  which  minimizes  the  overall  hot-  and  cold-clutter  power 


- 'H  zl 

wk  Rk  wk 


subject  to 


z'H  z l 

wk  Wfc 


1 


(24) 


is  given  by 

<o'k  = ^1(^+1)  (^r7)  (25) 

where  tti(-)  denotes  the  eigenvector  corresponding  to  the  minimum  eigenvalue  of  the  indicated  matrix.  We 
also  have  that  the  total  output  power  is  olnt  ~ <r^  + 

These  considerations  illustrate  two  important  issues.  Firstly,  when  the  specified  order  k is  sufficient  for 
effective  cold-clutter  mitigation  by  the  corresponding  rc-variate  MTI-type  filter,  “augmented”  MQ(rc+l) -variate 
STAP  can  simultaneously  reject  hot  and  cold  clutter.  The  hot-clutter  rejection  can  be  achieved  up  to  the  spatial 
(or  more  generally,  fast-time  spatio-temporal)  optimum  processing  limit,  while  the  cold-clutter  rejection  can 
be  achieved  up  to  the  temporal  optimum  processing  limit,  ie.  the  simultaneous  rejection  is  uncompromising  on 
its  effect  on  the  other  component.  Secondly,  each  of  the  («;+ 1)  blocks  of  this  filter  effectively  mitigates  the 
hot-clutter  component,  since  R ^ has  a block-diagonal  structure.  Indeed,  each  of  the  MQ-v ariate  components 
of  the  vector  wk  in  (25)  is  proportional  to  the  minimum  eigenvector  of  R^ . 

Now  suppose  that  the  total  sample  size  (T-Q+ 1)  used  for  (loaded)  sample  averaging  in  the  calculation  of 

- 2 

Rk  [2]  is  sufficiently  large  to  essentially  approach  the  efficiency  of  the  true  (deterministic)  value  for  the  STAP 
filter 

/zz  -t 

wfc  = (£fcJ  s{0 o)  (26) 


where  s(0 0)  = |sT(0o)  | O | . . . | oj  , then  we  expect  the  first  MQ -variate  section  of  this  “augmented” 

solution  to  provide  effective  fast-time  STAP  rejection  of  the  hot  clutter  for  the  kth  repetition  period,  it  was 
demonstrated  in  [2]  that  the  minimal  sample  volume  sufficient  for  3dB  average  losses  compared  with  the  ideal 
solution  (25)  double  the  signal  subspace  dimension  (ie.  rank)  of  the  noise-free  covariance  matrix.  Since  the 
covariance  matrix  rank  is  comparatively  small,  we  may  introduce  the  following  operational  routine. 


Algorithm  2 

Step  1 Given  the  observed  M-variate  snapshots  zkt  (k  = I ,.,N;  t = 1, . . . , T),  form  the  “stacked”  MQ-  ~ 
variate  samples  zkt  (k  = 1, . . . , AT;  t — 1, . . .,T  — Q+ 1),  and  the  “doubly  stacked”  MQ(k+1  ) -variate 
samples  zkt  ( k = k+1  t = 1, . . T-Q+l)  using(14). 

Step  2 Compute  the  (loaded)  sample  covariance  matrix  estimates 

1 T~Q+1  .K  zkH 
Rfc  jr ^ y ^kt  ® 

(27) 

for  k — k+  1, . . . , 2k+  1 and  some  small  constant  a,  then  compute  the  averaged  matrix 

zav  ^ K z % 

Rk+1  ~ , 1 ^K+j  ' 

7=0 

(28) 
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~ r - t 

Step  3 Let  Aq{0 o)  = A J ( 0O ) | O | . . . | O then  define  the  initial  (Ar  = /c+1)  stochastically  unconstrained 
and  range-independent  MQ(k+  l)-variate  fi  ter 


(cat/  \ “1  - ~ }f  / s av  \ — 1 _ 

R*+l  J Aq  (0O)  A?  (0O)  (R  K+lJ  Aq  W e?. 


The  first  M Q elements  of  this  vector  form  the  first  fast-time  STAP  filter  ie. 

v av  nT  = «t*  zo/w  I 

w*+l.t  (3°) 

where  = J/jwq  | O | . . . | oj  . 

Step  4 For  the  next  fast-time  STAP  filter  v\+2,t  (now  range-dependent  and  stochastically  constrained),  com- 
pute the  MQ(t c+  l)-variate  filter  w™+2  t by  applying  “sliding  window”  averaging: 

r*+2  = ^~j~Y  e R*+i+i  • (3i) 

j= o 

in  order  to  ensure  hot-clutter  rejection  over  the  repetition  periods  k = 2, . . . , /c+2. 

Compute  all  subsequent  operational  fast-time  STAP  filters  w™t  ( k = /c+2, . . . , N)  by 


- F) 


z av  'r  z av 

™k,  = ElMQwkt 

^y1  \-.zH  / s <w\  — i . , 1 — 1 

ftfc  ) -Afcf  /Ifcf  ( •Rjfe  ) -Ajbt  eJ+K 


where 


im  i * iz  -.j  r-^r  + 

**  = K |o|-|o] 


and  the  system  of  stochastic  constraints  is 


•~z  f . s(0o)  <\  - 

Akt  = Aq(0 o)  /mq  - Zkt  (35) 

where  Zfct  is  defined  by  (10). 


In  Step  3,  while  the  entire  MQ  (/c+1) -variate  filter  w™+1 1 would  provide  efifective  cancellation  for  both  hot 
and  cold  clutter,  we  expect  that  constructing  a fast-time  STAP  vector  from  its  first  MQ  elements  will 

deliver  equally  effective  cancellation  of  the  hot  clutter  only,  over  the  first  (/c+1)  repetition  periods.  This  will  be 
true  provided  that  the  number  of  degrees  of  freedom  in  the  vector  exceeds  the  rank  of  the  averaged  hot-clutter 

x av  w 

covariance  matrix  RK+1,  ie. 

MQ  > (/c+1)  [f(L+Q-1)]  (36) 

or  more  precisely,  taking  into  account  the  number  of  deterministic  and  stochastic  constraints: 


MQ  > (/c+1)  |p(Z,+Q-1)]  + (g+/c) . 


In  Step  4,  for  our  standard  example  of  k = 2,  once  again  the  initial  filter  w3t  rejects  hot  clutter  over 


■ 
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the  repetition  periods  k = 1,2,3,  while  the  second  filter  w™f  rejects  over  k — 2,3,4,  thus  ensuring  that  k 

2 Z 2 OH 

covariances  R.k  are  common  for  each  successive  average  Rk  . The  system  of  k stochastic  constraints  which 
ensures  the  stationarity  of  the  cold-clutter  component  is 

v>ZH  EMQ*kt  = zkt  = W™H  Emq  ikt  = v>7tH  Zkt  for  k = 2,  ...,k+1  (38) 

and  in  general 

•2  0.1 iH  _ = avH  _ c . 

wkt  zk-j,t  = V>kt  Zk—j,t  for  3 = (39) 

Thus  both  sides  of  these  constraints  will  consist  mainly  of  cold  clutter,  since  the  hot-clutter  component  is 
rejected. 

As  usual,  the  stochastic-constraints  approach  forms  the  kernel  of  this  algorithm.  The  main  difference  with 
all  previous  routines  is  the  method  of  separating  hot  and  cold  clutter:  here  we  embed  the  fast-time  STAP 
solution  into  the  more  general  “augmented”  3-D  STAP  one  to  extract  the  desired  solution  as  the  component  of 
the  MQ(k+ l)-variate  filter. 

It  should  be  clear  that  the  augmented  STAP  is  only  suitable  for  extracting  the  fast-time  STAP  solution,  and 
cannot  be  used  by  itself  as  a final  operational  routine  for  target  detection.  As  before,  this  is  because  the  very 
small  number  of  slow-time  lags  (k)  involved  in  cold-clutter  mitigation  leads  to  a wide  range  of  “blind”  Doppler 
frequencies,  where  any  targets  would  be  rejected  along  with  the  cold  clutter. 

3.  Operational  SC  STAP  Algorithm:  Simulation  and  Real  Data  Processing 

Naturally  both  of  the  proposed  operational  algorithms  need  to  be  verified.  The  most  important  question  for 
Algorithm  1 is  the  relevance  of  the  low-order  AR  cold-clutter  model,  and  this  shall  now  be  addressed  by 
real  data  processing,  since  there  is  no  point  filtering  simulated  cold  clutter  of  a known  AR  model.  The  other 
important  question  on  the  convergence  rate  of  the  process  to  the  true  value  as  T — >■  oo  for  Algorithm  1 has 
been  partly  addressed  in  our  previous  papers  on  supervised  training  [9,  1 8,  2]  with  the  established  sample  size 
T ~ 2(/v+ 1 )[P(L+Q  — 1)]  necessary  to  guarantee  average  losses  of  3dB  compared  with  the  optimal  solution. 

Thus  in  order  to  justify  the  first  operational  routine  we  have  to  demonstrate  that  via  a low-order  AR  model 
we  can  construct  a preprocessing  MTI-type  filter  that  can  efficiently  reject  cold  clutter. 'Firstly,  we  make  use 
of  the  same  real  SW  OTHR  data  as  in  [2],  where  the  CPI  is  100  seconds.  The  solid  line  in  Fig.  1(a)  shows  the 
standard  Doppler  spectrum  for  one  particular  range  cell;  note  that  the  subclutter  visibility  (the  main-peak-to- 
sidelobe  ratio)  is  approximately  50dB  in  this  case.  Forward  and  backward  averaging  has  been  used  to  define 
the  3x3(k  = 2)  inter-sweep  temporal  covariance  matrix  The  preprocessing  filter  has  been  defined  as 

= (£3)  e3  jef  (ilTj  e3|  (40) 

in  order  to  keep  unchanged  the  white-noise  output  power  (since  ||wa||  = 1).  The  dotted  line  in  Fig.  1(a) 
illustrates  the  Doppler  spectrum  of  the  residues  after  preprocessing.  We  see  that  the  “noise  floor”  obtained  by 
this  pre-filtering  is  essentially  the  same  as  in  the  initial  spectrum.  Moreover,  the  eigenvalues  of  the  sample 
matrix  7?-3 

Ai  = 3.42,  A2  = 0.0833,  A3  = 0.0020  (41) 

suggest  that  the  cold  clutter  could  be  rejected  by  about  27dB  compared  to  the  input  level.  This  agrees  with  the 
subclutter  visibility,  taking  into  account  the  compression  gain  of  N = 1000  repetition  periods  (~  30dB). 

Note  that  Fig.  1(a)  also  confirms  our  expectation  of  the  very  wide  “blind”  Doppler  bandwidth  which  makes 
this  type  of  processing  inappropriate  for  target  detection. 

Fig.  1(b)  illustrates  similar  processing  results  of  data  obtained  from  DSTO’s  HF  SWR  facility  located  at 
Port  Wakefield,  South  Australia.  This  radar  employs  a 16-element  linear  receiving  array  and  operates  over  the 
frequency  band  5 - 17MHz.  The  LFMCW  waveform  sweep  rate  and  bandwidth  are  selectable  over  a wide 
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range;  typically  the  waveform  repetition  frequency  is  around  4Hz  and  the  bandwidth  50kHz.  Since  the  repeti- 
tion period  is  almost  three  times  longer  than  in  the  previous  data,  the  energetic  components  of  the  sea-clutter 
Doppler  spectmm  (solid  line)  obviously  occupy  a significantly  wider  range  of  relative  Doppler  frequencies. 
Consequently,  the  third-order  preprocessing  filter  (k  = 2)  shown  by  the  dotted  line  is  not  extremely  effective, 
rejecting  the  most  prominent  Bragg  lines  barely  to  the  noise  floor  level.  The  “peak-to-noise”  ratio  at  the  output 
of  this  preprocessing  filter  (followed  by  a standard  weighted  FFT)  is  about  30dB,  while  the  initial  subclutter 
visibility  is  about  65dB.  Nevertheless,  if  we  increase  the  order  of  the  AR  model  to  k = 3,  then  the  correspond- 
ing fourth-order  preprocessing  filter  (dashed  line)  rejects  all  input  energetic  cold-clutter  components  far  below 
the  input  white  noise  level. 

Thus  for  typical  SW  OTHR  data,  small-order  AR  models  are  proven  to  provide  quite  effective  cold-clutter 
suppression,  in  turn  enabling  effective  hot-clutter-only  sample  extraction. 

Verification  of  the  second  operational  routine  is  not  so  simple,  but  the  following  two  important  questions 
may  be  addressed  by  simulation  studies.  Firstly,  wc  need  to  demonstrate  that  for  some  standard  hot-  and  cold- 
clutter  models,  hot-clutter  alone  can  be  rejected  by  an  MQ-\ ariate  component  of  the  “augmented”  MQ(k+ 
l)-variate  optimal  solution.  This  optimal  solution  is  constructed  from  the  exact  covariance  matrix,  and  is 
interpreted  as  the  limit  solution,  when  the  number  of  training  samples  tends  to  infinity.  Secondly,  we  need  to 
demonstrate  that  the  convergence  rate  with  respect  to  T is  sufficiently  high,  since  the  number  of  range  bins 
available  in  most  HF  OTHR  applications  is  usually  limited. 

The  following  simulations  arc  based  on  the  simple  scenario  of  pure  SAP  with  an  M — 16  element  antenna 
array  and  a single  fluctuating  jamming  source.  The  generalized  Watterson  model  described  in  [6,  2]  has  been 
used  to  simulate  the  spatial  and  temporal  (Doppler)  fluctuations  of  the  jammer.  The  spatial  and  temporal 
correlation  coefficients  have  been  chosen  to  reflect  typical  spatial  fluctuations  which  make  traditional  mitigation 
techniques  ineffective;  in  the  notation  of  [2],  these  parameters  are  Clf  — 0.90  and  pu  — 0.89  respectively,  with 
50dB  HCNR  and  9ie  = 59.1  degrees.  The  second-order  AR  model  with  N = 256  sweeps  per  dwell  and  the 
parameters  of  (11)  have  again  been  used  to  simulate  HF  scattering  from  the  sea  surface.  Fig.  2(a)  shows  the 
Doppler  spectrum  for  the  cold  clutter,  uncorrupted  by  any  jamming  signal,  at  the  output  of  the  conventional 
beamformer  (curve  labeled  “CBF  Y”).  We  see  that  the  subclutter  visibility  benchmark  is  approximately  80dB, 
which  is  an  upper  bound  for  most  practical  situations.  Also  illustrated  is  the  Doppler  spectrum  of  the  cold- 
clutter  signal  at  the  output  of  the  standard  (unconstrained)  optimal  SAP/STAP  filter: 

WST.AP  = 1 s(0 0)  (42) 

(note  that  here  R = R3^  and  s{6o)  — s(#o),  since  Q = 1)  (the  look  direction  is  0o  = 0).  Clearly  the  spatial 
nonstationarity  of  the  jammer  leads  to  a significant  degradation  in  subclutter  visibility  (~  -30dB),  due  to  the 
fluctuations  of  the  spatial  (antenna  array)  weight  vector.  For  this  simulation,  the  instantaneous  (per  sweep)  hot- 
clutter-to-cold-clutter  ratio  (HCCCR)  per  antenna  array  element  was  chosen  to  be  OdB,  so  that  the  subclutter 
visibility  for  the  conventional  beamformer  (“CBF  Z”)  is  also  about  50dB.  Therefore  in  this  particular  case, 
as  far  as  the  final  subclutter  visibility  is  concerned,  conventional  hot-clutter  mitigation  is  as  ineffective  as  no 
hot-clutter  mitigation  at  all. 

Fig.  2(b)  shows  the  results  of  our  analysis  into  the  limit  efficiency  of  the  algorithm  described  by  Eqns.  (32)— 

~Z  z.Z 

(35),  calculated  by  replacing  the  sample  matrix  Rk  by  the  true  covariance  matrix  Rk  , ie.  the  sample  volume  is 
infinitely  large.  The  curve  labeled  “long  Z”  presents  the  power  of  the  total  output  signal  for  the  “augmented” 
MQ(k+  l)-variate  STAP  filter  wkt\  “long  Y”  shows  the  power  of  the  cold-clutter  component  at  the  output  of 
this  filter;  “long  X”  corresponds  to  the  hot  clutter  power  for  the  same  filter,  while  “short  X”  illustrates  the  hot 
clutter  output  power  from  the  truncated  16-element  SAP  filter  wkt,  simply  formed  from  the  first  16  elements 
of  the  48 -variate  vector  wkt . 

We  see  that  in  the  given  scenario,  the  hot-clutter  component  is  rejected  far  below  the  level  of  the  cold- 
clutter  residues  at  the  output  of  the  augmented  STAP  filter  wkt.  Not  surprisingly,  the  short  filter  wk;  rejects 
the  hot  clutter  slightly  better  than  the  augmented  filter  due  to  the  independence  of  the  hot  clutter  over  adjacent 
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repetition  periods.  Thus  the  potential  effectiveness  of  the  proposed  routine  is  extremely  high,  since  the  scalar 
output  of  the  operational  filter  consists  almost  entirely  of  cold  clutter  only. 

The  most  important  remaining  question  deals  with  finite  sample  size.  In  this  regard,  Fig.  2(c)  differs  from 
the  previous  figure  only  by  the  use  of  the  finite  sample  size  T = 41.  For  this  reason,  the  “gap”  between  the  hot- 
and  cold-clutter  residues  is  smaller  than  in  the  previous  case  using  deterministic  covariance  matrix  calculations. 
Nevertheless,  this  gap  is  still  large  enough  to  guarantee  that  the  hot-clutter  component  at  the  output  of  the  filter 
wkt  will  be  negligible. 

Correspondingly,  the  Doppler  spectrum  of  the  total  signal  at  the  output  of  the  operational  filter  wkt  (curve 
labeled  “Op  SC  STAP”  in  Fig.  2(d))  is  practically  indistinguishable  from  the  cold-clutter-only  component 
Doppler  spectrum  at  the  output  of  the  conventional  bcamformcr  (“CBF  Y”),  and  from  the  spectrum  of  the 
signal  at  the  output  of  the  SC  STAP  filter  computed  with  the  true  covariance  matrix  (“Det  SC  STAP”). 

Thus  the  ability  of  the  operational  routine  to  suppress  fluctuating  jamming  signals  and  to  retain  the  initial 
subclutter  visibility  (very  high  in  this  example)  is  verified  in  this  case. 

In  order  to  explore  the  limitations  of  the  proposed  technique  within  the  framework  of  the  adopted  model,  we 
conducted  similar  simulations  for  input  HCCCRs  equal  to  1 OdB,  20dB  and  30dB  (Figs.  3, 4 and  5 respectively). 
We  see  that  only  in  the  (worst)  last  case  does  the  hot-clutter  output  power  approach  the  power  of  the  cold-clutter 
residues  at  the  output  of  the  augmented  adaptive  filter  w kt . Meanwhile,  subclutter  visibility  at  the  output  of  the 
operational  filter  wkt  is  degraded  up  to  the  level  of  about  70dB,  compared  with  the  initial  value  of  about  80dB. 

It  is  worth  mentioning  that  even  for  this  worst  case  (HCCCR=30dB),  the  potential  efficiency  of  hot-clutter 
mitigation  is  still  extremely  high.  This  means  that  with  an  appropriate  sample  volume,  one  can  approach  the 
initial  subclutter  visibility  level. 

Finally,  we  present  some  results  of  recent  field  trials  to  illustrate  the  efficiency  of  the  supervised  and  un- 
supervised SC  methods.  The  experimental  facility  involved  has  recently  been  developed  near  Darwin  by  the 
Australian  Defence  Science  and  Technology  Organisation  (DSTO),  in  collaboration  with  Telstra  Applied  Tech- 
nologies (TAT)  and  the  Cooperative  Research  Centre  for  Sensor  Signal  and  Information  Processing  (CSSIP); 
the  system  was  designed  to  function  primarily  as  a SW  OTHR,  though  it  can  operate  in  a variety  of  other  data 
acquisition  modes.  The  facility  consists  of  two  transmit  sites,  at  Stingray  Head  (65km  south-west  of  Darwin) 
and  Lee  Point  ( 1 Okm  north-east  of  Darwin),  together  with  a receive  site  at  Gunn  Point  (30km  north-east  of  Dar- 
win) [19].  The  receiving  system  is  based  on  a 32-element  uniform  linear  array,  some  500m  long,  connected  to 
a high  dynamic  range  32-channel  HF  receiver.  The  antenna  element  noise  figures  for  the  operational  frequency 
range  of  5-10MHz  do  not  exceed  5dB,  with  a figure  of  near  1 ,5dB  at  the  middle  of  the  band.  Such  an  antenna 
design  permits  a maximum  possible  extemal-to-intemal  noise  ratio,  and  thus  enables  efficient  external  noise 
mitigation. 

Our  results  illustrate  an  example  of  low-power  transmitter  operation,  where  the  last  30  or  so  range  cells 
of  the  80  available  ranges  have  the  sea-clutter  signal  deeply  submerged  into  the  environmental  noise.  There- 
fore we  were  able  to  use  these  last  30  ranges, for  supervised  training,  and  compare  the  efficiency  of  supervised 
and  unsupervised  training  SC  techniques  against  the  conventional  beamformer  (CBF)  and  standard  SAP  beam- 
former  averaged  over  the  entire  dwell.  For  unsupervised  training,  we  use  the  last  75  ranges,  since  the  first  five 
ranges  are  affected  by  the  extremely  strong  direct-wave  propagation. 

Algorithm  1 has  been  employed  with  the  30  range  cells  involved  in  cold-clutter  AR  parameter  estimation 
for  the  MTI-filter  design,  averaging  across  all  antenna  array  elements.  Figs.  6 — 8 present  the  detailed  range 
profiles  for  the  target  range  cell  19  and  37,  and  the  target-free  range  12.  The  total  noise  power  across  the 
frequency  bins  32  to  224  has  been  calculated  to  compare  the  signal-to-noise  ratio  improvement  (SNRI)  with 
respect  to  the  CBF,  bearing  in  mind  that  in  the  look  direction  (ideal  planar  wavefront),  all  beamformers  are 
normalised  to  the  same  gain.  Interestingly,  this  analysis  demonstrates  that  unsupervised  training  (10-1 5dB 
SNRI)  is  slightly  better  than  supervised  training  (8-13dB  SNRI)  in  this  case.  This  could  be  explained  by  an 
improved  interference  averaging  with  most  of  the  available  ranges  involved  (75  out  of  80).  Practically,  though, 
this  improvement  is  not  indicated  by  the  actual  SNRIs  for  the  reference  targets.  This  discrepancy  between  the 
expected  and  actual  SNRIs  is  once  again  explained  by  the  fact  that  the  MT1  residues  that  contain  the  target 
components  have  been  used  for  interference  covariance  matrix  estimation.  Since  perfect  antenna  calibration  is 
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not  attainable  in  practice,  some  degradation  in  SNR1  is  inevitable.  This  is  a well-known  phenomenon,  and  a 
straight-forward  modification  that  excludes  target-suspicious  range  cells  from  averaging  provides  one  antidote. 

Nevertheless,  both  SC  SAP  options  demonstrate  the  considerable  improvement  achievable  in  practice  com- 
pared with  CBF  and  averaged  SAP  in  this  environmental  noise  situation. 


4.  Summary  and  Conclusions 

We  have  extended  the  domain  of  practical  application  of  the  stochastic-constraints  (SC)  method  to  unsupervised 
training  scenarios,  typical  of  existing  FMCW  OTHR  systems. 

Two  operational  routines  have  been  proposed  here.  The  first  technique  includes  “slow-time”  preprocessing 
of  the  input  data  by  the  MTI-type  filter  that  can  suppress  the  cold-clutter  component  far  below  the  jamming 
signal  level.  We  have  demonstrated  by  real  SW  OTHR  data  processing  that  effective  cold-clutter  rejection  can 
be  achieved  by  involving  a very  modest  number  of  repetition  periods  in  preprocessing.  In  fact,  the  order  of  the 
preprocessing  filter  is  equal  to  the  order  of  the  AR  cold-clutter  model,  which  in  turn  is  equal  to  the  number 
of  stochastic  constraints  that  secure  the  stationarity  of  the  output  cold-clutter  signal.  It  is  highly  significant 
that  low-order  AR  models  have  now  been  proven  to  be  adequate  only  for  the  local  description  of  the  cold 
clutter,  ie.  over  a very  small  number  of  consecutive  repetition  periods.  Correspondingly,  the  preprocessing 
approach  that  we  have  introduced  can  only  be  used  for  the  extraction  of  the  hot-clutter  signal,  and  is  completely 
inappropriate  for  taiget  detection  due  to  the  unacceptably  broad  “blind”  Doppler  frequency  bandwidth. 

Within  this  first  approach,  the  AR  model,  and  consequently  the  preprocessing  filter,  are  assumed  to  be 
known  a priori  (or  estimated).  The  second,  more  general  approach,  relies  only  upon  the  chosen  order  of  pre- 
processing filter  sufficient  for  cold-clutter  rejection,  while  the  second-order  moments  for  both  hot  and  cold 
clutter  are  unknown.  In  this  second  method,  the  desired  fast-time  SAP/STAP  hot-clutter  rejection  filter  is  de- 
fined as  part  of  the  “augmented”  STAP  filter  which  involves  (k+1)  consecutive  repetition  periods  for  additional 
cold-clutter  mitigation.  More  precisely,  the  MCJ-variate  fast-time  operational  STAP  filter  (or  Af-v ariate  SAP 
filter)  Wfr  is  defined  as  the  vector  consisting  of  the  first  MQ  (or  M)  elements  of  the  augmented  MQ(k+  1)- 
variate  STAP  filter  (or  the  M(k+ 1) -variate  STAP  filter)  w^. 

While  this  augmented  filter  ibkt  simultaneously  provides  both  hot-  and  cold-clutter  rejection,  the  short 
version  wkt  was  demonstrated  to  reject  effectively  only  the  hot-clutter  component,  far  below  the  cold-clutter 
signal  level.  The  proposed  stochastic-constraints  method  ensures  the  stationarity  of  the  scalar  cold-clutter  signal 
at  the  output  of  the  operational  filter  wkt . 

Simulation  results  presented  demonstrate  that,  for  a typical  scenario,  the  SC  technique  provides  effective 
jammer  mitigation,  retaining  the  initial  high  level  of  subclutter  visibility.  These  simulations  involved  a quite 
modest  number  of  range  resolution  cells,  typical  of  existing  FMCW  OTHR  installations.  The  boundaries  of 
usefulness  of  the  SC  method  are  explored  in  terms  of  subclutter  visibility  as  a function  of  input  hot-  to  cold- 
clutter  ratio.  Most  of  the  efficiency  analysis  results  presented  here  are  obtained  by  direct  simulations,  while  a 
proper  analytic  study  of  the  convergence  properties  of  the  SC  approach  remains  to  be  done.  Nevertheless,  we 
have  demonstrated  by  real  SW  OTHR  data  processing  that  unsupervised  training  scenarios  typical  of  FMCW 
OTHR  can  be  successfully  treated  in  an  operational  mode  using  stochastic-constraints  principles.  - 

Naturally  the  techniques  introduced  here  could  be  useful  for  applications  other  than  HF  OTHR.  The  simi- 
larities between  HF  OTHR  and  airborne  radar  have  already  been  discussed  in  [2]. 
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Figure  1 : Noise-floor  comparison  in  two  sets  of  real  surface-wave  OTHR  data. 


Figure  2:  Simulation  results  for  a hot-clutter-to-cold-clutter  ratio  (HCCCR)  of  OdB:  (a)  standard  processing 
Doppler  spectra,  (b)  and  (c)  scalar  filter  output,  (d)  Doppler  spectra  for  proposed  techniques. 


= 0.89,  £ = 0.9,  HCCCR  = 20dB  M = 16,  N = 256,  T = 4I.  Q = l,  p = 0.89,  £ = 0.9,  HCCCR  = 20dB,  deterministic  R 


repetition  period  (d)  relative  Doppler  frequency 


0.89,  £-0.9,  HCCCR-30dB  M-16,  N-256,  T - 41 , Q-l,  p-0.89,  5*0.9.  HCCCR  - 30dB,  deterministic  R 
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(c)  repetition  period  (d)  relative  Doppler  frequency 


